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■ Abstract. We present rigorous estimates for some physical quantities related to 
turbulent and non-turbulent channel flows driven by a uniform pressure gradient. 

• Such results are based on the concept of stationary statistical solution, which is 

Q \ related to the notion of ensemble average for flows in statistical equilibrium. We 

provide a lower bound estimate for the mean skin friction coefficient and improve 
on a previous upper bound estimate for the same quantity; both estimates in terms 
of the Reynolds number. We also present lower and upper bound estimates for the 
^—1 \ mean rate of energy dissipation, the mean longitudinal velocity (in the direction of 

J> ■ the pressure gradient), and the mean kinetic energy. In particular, we obtain an 

upper bound related to the energy dissipation law, namely that the mean rate of 
energy dissipation is essentially bounded by a non-dimensional universal constant 
times the cube of the mean longitudinal velocity over a characteristic macro-scale 
^—1 ■ length. Finally, we investigate the scale-by-scale energy injection due to the pressure 

\Q \ gradient, proving an upper bound estimate for the decrease of this energy injection 

■ as the scale length decreases. 

1. Introduction 

CZ • 

While the existence of exact solutions of the Navier-Stokes Equations are not avail- 
J> ! able in general, most of the classical research on turbulence theory consist of ap- 
proximate methods based on a few exact deductions, supplemented with intuitive 
hypotheses about the nature of the phenomenon, such as scaling assumptions and 



moment truncation models; see for example [H ■ 

Recently, part of the theoretical research on turbulence has been concentrated 
on deriving rigorous bounds on characteristic quantities of turbulent flows directly 
from the equations of motion. These results are important to substantiate the ones 
obtained via the classical approximation methods. 

Decomposing the turbulent flow into a stationary background flow and a fluctu- 
ation component, and using variational methods, Constantin and Doering derived 
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rigorous results for the long-time averaged rate of energy dissipation of flows in some 
geometries, in particular for the channel flow driven by a pressure gradient, in which 
case the estimate also yields an estimate for the mean skin friction coefficient; see 



Meanwhile, rigorous results were recently established for the three-dimensional the- 
ory of homogeneous stationary statistical turbulence in [131 [HI [15] using the concepts 
of stationary statistical solutions of the Navier-Stokes equations and generalized time 
average measures, and using energy-type methods. 

This paper presents a combination of those results in the specific case of channel 
flows driven by a uniform pressure gradient. More specifically, we extend the upper 
bound estimate for the long-time averaged rate skin friction coefficient, obtained in [4], 
to general stationary statistical solutions, slightly simplifying their proof and slightly 
improving their estimates. We also obtain a lower bound estimate for the skin friction 
coefficient, which cannot be obtained by the variational principle method of [4]. More 
precisely, we show that for every stationary statistical solution the corresponding skin 
friction coefficient Cf satisfies 



for low Reynolds number flows, and 



for high Reynolds number flows, where the Reynolds number is defined by Re = hi// u, 
with h being the height of the channel and U, the mean longitudinal velocity. 

The lower-bound estimate for Cf is nearly optimal in the sense that the stationary 
statistical solution is arbitrary and may be concentrated on the plane Poiseuille flow 
(which is unstable for high Reynolds number flows, but anyway exists in a mathemat- 
ical sense), for which Cf = 12/Re. The upper-bound estimate might not be optimal 
since heuristic arguments suggest that Cf ~ (In Re) -2 for high- Reynolds number tur- 
bulent flows. Nevertheless, it represents a nearly 19% improvement over the estimate 
obtained in [4] on the leading order constant term (from 0.597 to 0.484.) 

We also give upper and lower bound estimates for some other physical quantities, 
such as the mean energy dissipation rate, the mean kinetic energy, and the mean 
longitudinal velocity. In particular, we prove an upper bound estimate related to the 
energy dissipation law, namely that for large Reynolds number flows the mean rate 
of energy dissipation is essentially bounded by a non-dimensional universal constant 
times the cube of the mean longitudinal velocity U over the height h of the channel: 



The leading order constant term obtained in [4] was approximately 0.0884, and it was 
remarked in that work that this term is much lower than 1, hinting that this result 



13. 4j. 
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is substantially more than a formalized dimensional analysis argument. The same 
applies here. 

Finally, we study the scale-by-scale energy injection term due to the pressure gra- 
dient. We show that the energy injected into the modes larger than or equal to k 
is bounded by a term proportional to k~ 3 / 2 . The motivation for the study of the 
decrease of energy injection comes from the Kolmogorov theory of turbulence. This 
theory argues that for turbulent flows there is a certain range of scales much lower 
than the energy injection scales and greater than the energy dissipative scales in which 
the kinetic energy is transferred to the small scales at a nearly constant rate equal 
to the energy dissipation rate. This theory was proposed in the idealized case of lo- 
cally homogeneous turbulence, away from the boundaries, under the assumption that 
the energy injection is concentrated on the large scales. However, it is known from 
experiments that for wall bounded turbulence this hypothesis needs to be corrected 
[H [281 [23]. In particular, the energy injection occurs at arbitrarily small scales. The 
estimates presented yield an upper bound on the rate of decrease of energy injection 
as the scale length decreases. 

The remaining of the paper is organized as follows. In the next section we introduce 
the convenient mathematical setting used throughout the work. In Section 3, we 
define the notion of stationary statistical solution of the Navier-Stokes equations and 
present some related results. In Section 4, we rigorously define the characteristic 
quantities that will be estimated, such as the mean energy dissipation rate, mean 
kinetic energy, mean longitudinal velocity, and mean skin friction coefficient. In 
Section 5, we establish a relation between stationary statistical solutions and time 
averages. In Sections 6 and 7, we explicitly derive rigorous bounds for the mentioned 
physical quantities, utilizing both methods of [131 [HI [15] an< ^ °^ @J- m Section 8, we 
list the explicit values of the mentioned physical quantities for the specific case of the 
laminar Poiseuille flow, verifying that some of the results obtained in Section 6 and 
Section 7 are optimal in some sense. In Section 9, we conclude the work with the 
discussion about the scale-by- scale energy injection. 

2. Mathematical framework of the Navier-Stokes Equations 

We consider an incompressible Newtonian flow confined to a rectangular periodic 
channel and driven by a uniform pressure gradient. More precisely, the velocity vector 
field u = (ui,U2, U3) of the fluid satisfies the incompressible Navier-Stokes equations 



in the domain Q = (0, L x ) x (0, L y ) x (0, h). The scalar p is the kinematic pressure. 
We denote by x = (x, y, z) the space variable. The boundary conditions are no-slip 
on the planes z = and z = h and periodic in the x and y directions, with periods L x 
and L y , respectively, for both u and p. The parameter P/L x denotes the magnitude 
of the applied pressure gradient. The parameter v > is the kinematic viscosity, e 1 




V ■ u = 0, 
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is the unit vector in the x direction, and L x , L y , h, P > 0. We sometimes refer to the 
direction x of the pressure gradient as the longitudinal direction. 

The mathematical formulation of the Navier-Stokes equations in this geometry can 
be easily adapted from the no-slip or fully-periodic case developed in [5| \13 \ [T9l 1251 [27] . 

The formulation yields a functional equation for the time-dependent velocity field 
u = u(t) of the form: 



du . . 

dT = F < u > 



fp — vAxx — B(u, u) 



where 



(2.2) 



(2.3) 



Two fundamental spaces are defined by 



H 



wG(lL(R 2 x(0,/i))f, V-w = 0, 



u 



w n ; 



w(x + L x ,y, z) 
w(x,y + L y ,z) 



w(x,y,z), 

w(x,y,z), a.e. (x,y,z) G 



x (0,h). 



> . 



w 3 (x, y, 0) = w 3 {x, y, h) = 0, a.e. (x, y) G 



and 



V 



u 



wG(iJi c (M 2 x(0,/i))) 3 , V-w = 0, 
w(x + L x ,y, z) = w(x,y,z), 
w(x,y + Ly,z)=w(x,y,z), a.e. (x,y,z)e 
w(x, y, 0) = w(x, y, h) = 0, a.e. (x, y) G M 2 

The inner products in H and V are denoted respectively by 

du 



x (0,h). 



u, V 



u(x)-v(x)dx, ((u,v) 



1 1/2 



E 

n t=l,3 



vl/2 



(9v 



dx, 



and the associated norms by |u| = (u, uj"'", ||u|| = ^u, uj y 

We identify if with its dual and consider the dual space V of V, so that V C if C 
V, with the injections being continuous, each space dense in the following one. We 
also denote by H w the space H endowed with its weak topology. 

We denote by Plh the (Leray-Helmhotz) orthogonal projector in L 2 (Q) 3 onto the 
subspace H. The operator A in (I2.2p is the Stokes operator given by Au = — P L h^ u - 
The term B(u, v) = Plh((u • V)v) is a bilinear term associated with the inertial 
term. Moreover, since the Stokes operator is a positive self-adjoint operator on H, 
we consider its powers A s , s G R, with domain D(A S ). We have V = D(A 1/2 ) and its 
dual V = D{A^I 2 ). 
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The Stokes operator possesses a complete orthonormal basis of eigenvectors in H, 
{Wj,i,k}j,i,k, o f the form 

w jtl}k (x,y,z) = exp ^tt + J^j ™j,iA z ), ( 2 - 4 ) 

where (j,l,k) e Z x Z x N, -Aw^y^ = \j,i,k™j,i,k, and each w^^z) satisfies a one- 
dimensional eigenvalue problem, with < Xj t i t k — * oo, when j, I, k — > oo. We write 
the spectral expansion of u in this basis as 

u(x,y,z) = ^2uj,i, k ^j,i,k{x,y,z). (2.5) 

j,l,k 

To each eigenvalue Xj,i,k we associate a wavenumber k = /Cj,j,fc = A^ 2 fc . Since u e V 
vanishes on the top and bottom walls, Poincare inequality applies, yielding a bound 
on |u| in terms of ||u||. In fact, we have precisely 

|u|o < Ar 1 ||u|| 2 , (2.6) 
where Ai = 7i 2 /h 2 is the smallest positive eigenvalue of the Stokes operator on this 

1 /2 

geometry. The smallest positive wavenumber is K\ = 

Ar = n/h. 

We define the component u K of the vector field u, for a single wavenumber k, by 

u « = E v-j,i,kWj,i,k, 

and the component u K / jK » with a range of wave numbers [k 1 , k") by 

U K ' j/c " = E u K . 

k'<k<k" 

We then write the Navier-Stokes equations projected on those components in the 
form 

^f- + vAu K ,, K „ + B(u, n)^ = (fp) K , K „ , (2.7) 

where 

(fp) K = E (fp.w iiI>fc )w WiJfc . (2.8) 

Kj j,l,k =K 

Now, taking the inner product in H of the bilinear term with a third variable yields 
a trilinear term 

b(u, v, w) = (-B(u, v),w), 
which is defined for u, v, w in V. An important relation for the trilinear term is the 
orthogonality property 

6(u,v,v) = 0, (2.9) 
for u, v G V. It follows from this relation the anti-symmetry property 

6(u, v, w) = -6(u, w, v), (2.10) 

for u, v, w e V. 
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2.3 



3. Statistical solutions and the Reynolds Equations 



A mathematical framework for the conventional theory of turbulence is based on 
the concept of stationary statistical solution of the Navier-Stokes equations. This 
amounts to considering the space H as a probability space with the cx-algebra of 
the Borel sets of H and endowed with a Borel probability measure. The ensemble 
averages are then regarded as averages with respect to this Borel probability measure. 
In our three-dimensional case we work mostly with the weak topology. Fortunately 
the Borel <r-algebra generated by the weakly open sets coincides with that for the 
open sets in the strong topology. Since H is a separable Hilbert space every Borel 
measure is automatically regular. An important consequence of the regularity of a 
Borel probability measure is the density of the continuous functions (or just weakly 
continuous functions) in the space of integrable functions. 

We say that a measure \x in H is carried by a measurable set E when E has full 
measure in H, i.e. // (H \ E) = 0. The support of a Borel probability measure \i is the 
smallest closed set which carries /i. The ensemble averages are regarded as averages 
with respect to a Borel probability measure fi on H. If ip : H — > R is a Borel function 
representing some physical information (p (u) extracted from a velocity field u, such 
as kinetic energy, velocity, enstrophy, etc., then its mean value is 



The reader is referred to [13] for more details. 

Now, we define a class of Borel functions that are particularly useful in order to 
make a rigorous definition of a stationary statistical solution of the Navier-Stokes 
equations. 

Definition 3.1. We define the class T of test functions to be the set of real-valued 
functionals \I> = ^(u) on H that are bounded on bounded subsets of H and such that 
the following conditions hold: 

(1) For any u G V , the Frechet derivative ^/'(u) taken in H along V exists. More 
precisely, for each u G V , there exists an element in H denoted \&'(u) such 
that 



(2) ^/'(u) G V for all u G V , and u — > \I/'(u) is continuous and bounded as a 
function from V into V . 

For example, we can take the cylindrical test functions : H — > R of the form 
^(u) = ip ((u, gi), . . . , (u, g m )), where ip is a C 1 scalar function on R m , m G N, with 







STATISTICAL ESTIMATES FOR CHANNEL FLOWS DRIVEN BY A PRESSURE GRADIENT 7 



compact support, and gi, . . . , g m belong to V. For this case we have 

m 

= d ^^ u ' ■ ■ ■ , (u, g J ))g„ 
i=i 

where denotes the derivative of ip with respect to the j-th variable. It follows 
that \P'(u) G V since it is a linear combination of the g 3 -. 

Now, we define the notion of a stationary statistical solution of the Navier-Stokes 
equations. 

Definition 3.2. A stationary statistical solution of the Navier-Stokes equation is a 
Borel probability measure fi on H such that 

(1) / ||u|| 2 d/i(u) < oo; 
Jh 

(2) / (F(u), *'(u))d/i(u) = 0, for any *eT, where F(u) is as in QZSJ); 

(3) / {z/ ||u|| 2 — (fp, u)} <i/i(u) < 0, for all < ei < e 2 < +oo. 

iei<|u|2/2<e 2 

The first condition means that an arbitrary stationary statistical solution has finite 
mean enstrophy. This is natural when we compare with individual solutions whose 
time average is bounded uniformly with respect to the time interval. It is also needed 
to make sense out of the second condition. 

The last condition on the definition above is an energy-type inequality, and one 
can deduce from it that the support of a stationary statistical solution is included in 
the weak attractor A w , see [HUE]) which is bounded in H according to 

uG* uh 1 / 2 

\<< — = ^j^G\ VueX, (3.3) 

Cr — 9 : /9 \A Ip o, 
V TV ' 

where G* is a nondimensional number called Grashof number. 

The concept of stationary statistical solution is regarded as a generalization of the 
notion of invariant measure. It is relevant to our three-dimensional case, in which a 
semigroup is not well-defined. 

Due to these regularity properties obtained for stationary statistical solutions (finite 
mean enstrophy and with support bounded in H), the mean value (<£>(u)) can be 
defined not only for weakly continuous functions bounded in H but for any real- 
valued function if which is continuous in V and satisfies the estimate 

\<p\ <C(\u\ )(l + v- 2 \\ /2 \\u\\\ (3.4) 

where C(|u| Q ) is bounded on bounded subsets of H. Important examples of such ip 
are |u|J, ||u|| 2 , 6(u KOjK , u KOiK , u Kj00 ), and b(u K 
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By a duality argument we can extend the ensemble averages to functions with value 
in some function space. More precisely, we define the velocity field (u) and the mean 
value (B(u, u)) of the inertial term by 

«u),v)= / (u,v)^(u), WveV, 
Jh 

«B(u,u)),v)= / (5(u, u ),v)^(u), WED (A 3 / 8 ). 
Jh 

The mean flow (u) is a vector field on Cl with (u) £ V, while {B(u, u)) £ _D(v4~ 3 / 8 ). 

Since we assume statistical equilibrium, the stationary form of the Reynolds equa- 
tions can be recovered within this framework; see also 



Proposition 3.3. Given a stationary statistical solution in the sense of Definition 
\3.S\ the following functional form of the Reynolds equations hold in V: 

uA(u) + (B(u,u)} = f P . (3.5) 

Proof. Let ip be a C l real-valued function with compact support on R. For any 
v £ V and any wavenumber k, the function $(u) = if>((u, v K1 „)) is a cylindrical test 
function. Thus, 

/ ip'((u, v KlA )) {(f P , v KljK ) - u(Au, v K1)K ) - b(u, u, v KljK )} dfi(u) = 0. 
Jh 

Let if}' converge pointwise to 1 while being uniformly bounded, so that at the limit 
we find 

/ {(fp, v K i,k) - v{Au, v KuK ) - 6(u, u, v Kl . K )} d//(u) = 0. 
Jh 

For each fixed v £ V, we may let k go to infinity to find (since fi has finite enstrophy 
and with support bounded in H) 

I {(fp, v) — u(Au, v) — 6(u, u, v)} d/i(u) = 0. 
Jh 

which gives us the result. □ 

We end this section with a result concerning the Grashof number G*, which yields 
a bound in H on the weak attractor A m in terms of P: 



Lemma 3.1. We have, more explicitly, 

G* = W?» P. (3.6) 



Proof. Since fp = (P/L x )e-j_, we have that 



A- 1 f P =(^(/i-z),0,0). 
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Hence, 

\A-Wf P \l = (fp, A-%) = [ ^^z(h - z)dx = ^f P \ (3.7) 

Taking the square root of the equality above and substituting in the definition of the 
Grashof number give us the result. □ 

Remark 3.1. The vector-field A~ l ip is directly related to the plane Poiseuille flow. 
In fact, the plane Poiseuille flow is precisely u = A~ 1 fp/v = (Pz(h — z)/2uL x ,0,0); 
see Section 7. 

4. Characteristic dimensions and nondimensional numbers 

The macroscopic characteristic length is considered to be h and the macroscale 
characteristic wavenumber is hq = 1/h. The total mass of the fluid in the channel 
is p Q L x L y h, where po denotes the uniform mass density of the fluid. Then, for a 
given stationary statistical solution /i, the corresponding mean kinetic energy per 
unit mass and the mean energy dissipation rate per unit time and unit mass are 
given respectively by 

*2L x Lyh L x Lyli 
The mean longitudinal velocity is defined by 

U = Tli / (/ / u ^ x iVi z ) d y d ^j dfx(u). (4.1) 

Note that this definition makes sense and the expression does not depend on x due 
to the incompressibility and boundary conditions. 

With this velocity scale we may define the following Reynolds number 

Re = ™. (4.2) 
v 

A dimensionless ratio of the applied pressure gradient to the square of the flow velocity 
scale is provided by the skin friction coefficient given by 

<"-£*■ 

Thanks to condition (3) of Definition 13.21 and to the divergence-free and boundary 
conditions, the mean longitudinal velocity and the mean energy dissipation rate are 
related by 

« < U ~f- («) 

i-'x 

Now, suppose that u(x, t) is a weak solution of the Navier-Stokes equations (12.11) 
with initial condition u (x). We define the finite-time average longitudinal velocity 
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by 



1 1 



Ut — 7Tm / I / / Ui(x,y, z,t)dydz ) dt. (4.5) 



L y h T j \j j 

This expression is also well defined and independend of x due to the incompressibility 
and boundary conditions. 

In the sequel, we will also consider the similarly defined time-averaged dissipation 
rate 

v 9, v 1 



and time-averaged kinetic energy 

2L x L y h 2L x L y h T J 

Notice that even if finite-time averages are bounded, their long-time limits may not 
exist. 

5. Time averages and stationary statistical solutions 

Since the usual limit of long-time averaged quantities may not exist, we aim to 
obtain eventual bounds for these time averaged quantities. In this section, we will 
establish, via generalized limits, a rigorous relationship between certain limits of these 
quantities and the stationary statistical solutions, see [13] . 

For example, suppose we are interested in estimating the upper limit of the time 
averaged velocity of a weak solution u(x, t) of (12. ip 

l i r T ( r L * r h \ 

U = — - lim sup — / I / / ui(x,y, z,t)dydz \ dt. (5.1) 
Lyh t^oo T J \J Q J J 

Since an upper bound for the mean longitudinal velocity U, associated with an arbi- 
trary stationary statistical solution /x, is derived in (16. 1ft . namely 

U < ^P, 
1 QuL x 

we may establish a direct relation between the time average longitudinal velocity (15.11) 
and the mean longitudinal velocity associated with a specific stationary statistical 
solution /Jo, such as 

lim sup — [ ( I / ui{x,y, z,t)dydz\ dt = ( I Ui(x,y, z)dydz\ d/^o, 
t^oo T Jo \J Jo J J H \Jq Jo J 

(5.2) 

in such a way that we can give an upper bound to ( 15.11) using (16.11) . obtaining 

lim sup — J I J J u 1 (x,y,z,t)dydz\dt=U l t <-^-j—P. (5.3) 



Lyh 
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This relation between long-time averages and stationary statistical solutions is re- 
alized via the notion of generalized limit, which is defined as follows 

Definition 5.1. A generalized limit is any linear functional, denoted Lim^oo; de- 
fined on the space B([0, <x>)) of all bounded real-valued functions on [0, oo) and satis- 
fying 

(1) LiMt_>oo £(r) > 0, Vg G B([0,oo)) with g{s) > 0, Vs > 0; 

(2) LiMt-+oo g(T) = \hn T -. o0 g(T), Vg G B([0, oo)) 

such that the classical limit, denoted limj^oo, exists. 

Remark 5.1. It can be shown that given a particular g G B([0, oo)) and a sequence 
tj — > oo for which go converges to a number I, there exists a generalized limit LlMr^oo 
satisfying LiMt^oo #0 = I', see [2| [13]. 

Proposition 5.1. Let (p G C(H W ). Suppose that for every stationary statistical 
solution \i, the associated average of (p satisfies (<p) < C\, for some constant C\. 
Then, given a weak solution w(x, t) defined on [0, oo), we have 

i r T 

limsup-/ (p(w(t))dt < C x . (5.4) 

Similarly, if for some constant C 2 we have(ip) > C 2 for every stationary statistical 
solution, then 

liminfi f <p(w(t))dt > C 2 . (5.5) 

Proof. We will prove inequality (15 .4p . Inequality (15.51) follows by a similar argument. 
Let w = w(0). Consider the set 

K w = {v G H; |v| 2 < \w \l + \f P \ 2 /u 2 \\} , 

endowed with the weak topology of H . K w is compact in H w and is such that 
w(t) G K w , for all t > 0; see [SJ 

Let t/) G C(K w ). Since l"f w is compact, the function £ i— > ip(w(t)) is continuous and 
bounded. Thus, 

1 " 



0o(*) = ~ ip(w(s))ds 

£ JO 

makes sense, and is continuous and bounded for t > 0. Therefore, its generalized 
limit is well defined, and by Remark 15.11 if we choose a subsequence tj — > oo for 
which go(tj) converges to limsup^^ go(t), there exists a generalized limit LlM T 
satisfying 



1 /" T 1 /■* 

LlM — / -0(w(t))c?t = limsup — / 

T^oo T J t^oo t J 



rT 

if)(w(s))ds. 
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Now, we relate this generalized time average with stationary statistical solutions. 
Since the weak solution 1 1— > w(t) belongs to the compact set K w in H w , see [13], and 
since 

i r T 

i) i-> Lim - / V(w(t))d* 

is a positive linear functional on C(K W ), we use the Kakutani-Riesz representation 
theorem, see [29], and conclude that there exists a measure fi on H such that 

LlM^/ V(w(x,t))dt= / ip(u)dfi (u), (5.6) 

for all t/) G C(K w ). It is shown in [13] that /^o defined above is a stationary statistical 
solution. 

Therefore, since ip\ic w £ C (if™), and /i(H\K w ) = 0, for every stationary statistical 
solution, see [13], and in particular for /i , we conclude that 

limsup- / if(w(x, t))dt — I (p(u)du, (u). 

t^oo t Jq J h 

Thus, since (ip) < C\ for all stationary statistical solution /z, and in particular /i , we 
have 

limsup— / ip(w(t))dt = / ip(u)dfio(u) < C\. □ 
r^oo T J J h 

Now, returning to the mean longitudinal velocity example, since 

I / ui(x,y, z)dydz 



belongs to C(K W ), there exists a stationary statistical solution fi satisfying (15. 2L 
which together with (16. ip . yields the upper bound (15. 3ft . 

Remark 5.2. Proposition 15.11 shows that every estimate involving the average of a 
continuous quantity on C(K W ) can be stated as a superior or inferior limit of its time 
average. 

This is true for the energy injection term, (fp, u), and also for the mean longitudinal 
velocity U. However, we are also interested in estimating quantities involving |u| 
and ||u||, which are not weakly continuous. Fortunately, we are still able to estimate 
these quantities by approximating via Galerkin projections as shown in the next 
proposition. 

Proposition 5.2. Let w(x, t) be a weak solution of the NSE defined on [0,oo), and 
suppose that for every stationary statistical solution fi, we have the following bounds 



G x <\ |u|^/i (u) < C 2 , 
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and 



C 3 < ||u|| 2 d/i (u). 
Jh 



Then, we also have 



and 



liminf i / \w(t)\ 2 dt>C u (5.7) 

T^oo 1 J Q 

i r T 

lim sup- / \w{t)\ 2 dt < C 2 , (5.8) 

T^oo J- Jo 

i r T 

liminf- / \\w(t)\\ 2 dt > C 3 . (5.9) 



Proof. Since |-P K U I > II-Pk u II e C(K W ), where P K are the usual Galerkin projectors, we 
have by (15.6p that given a stationary statistical solution /i generated by a generalized 
time average, the following equations are valid 

-T 



and 

Now, since 
and 



LlM^/ \P lt w(t)\ 2 dt = [ \P K u\ 2 d/i (u), (5.10) 

T^oo ± J Q J H 

LlM^/ \\P K w(t)\\ 2 dt= [ \\P K u\\ 2 dfi (u). (5.11) 



\w(t% - |P,w(f)|£ = |Q„w(f)|£ < ^ 2 ||w(t)| 

fT 



w(t)|| 2 dt < C < oo, 



1 

T 

where C is independent of T [5J[2S], we have by the usual properties of the generalized 
limits that 

i pT i r T i r T 

Lim — / |w(t) | 2 (it — LlM — / \P K w(t)\ 2 dt < k~ 2 LlM — / || w(t) || 2 c?t 

1 r T 

<fir 2 limsup-/ \\w(t)\\ 2 dt < Ck~ 2 -> 0, k -> oo. (5.12) 

Thus, considering the generalized limit, LiMt^oo, that extends the left hand side of 
(15.71) . we have 

liminf i f \w(t)\ 2 dt= Lim! / |w(t)| 2 dt = lim Lim ± / |P K w(t)| 2 rft 



= lim / |P K u| dfMo(u) = / |u| g^o( u ) > Ci, (5.13) 
J// Jh 

where the last equality of the expression above follows from the Monotone Conver- 
gence Theorem. The bound (I5.8P follows in a similar way. 
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Now, we will prove ( I5.9p . Consider the generalized limit, LiMt^oo, that extends 
the l.h.s. of (15. 9p . and notice that 

||w(i)|| > \\PMt)\\- 

Then, we have 

liminfi / \\w(t)f Q dt= LiM^ / \\w(t)f dt> lim LlM ^ / ||P K w(t)|| 2 dt 

= lim / ||P K u|| 2 ^ (u) = / ||u|| 2 ^ (u) > C 3 , (5.14) 

where, again, the last equality in the expression above follows from the Monotone 
Convergence Theorem. □ 

Remark 5.3. Propositions 15.11 and 15.21 show that, except for Theorem 17.81 and 
Propositon I7.9[ every estimate in the sequel can be stated as a superior or inferior 
limit of their time averages. The reason why these results do not apply to Theorem 
17.81 and Propositon 17.91 is that they involve an upper bound to ||w||, which is not 
considered by the propositions above. However, they can still be stated in terms of 
their time averages as seen in Remark 17.11 



6. Estimates on the mean longitudinal velocity and on the skin 

friction coefficient 

We start by deriving an upper bound on the mean longitudinal velocity. 

Theorem 6.1. For every stationary statistical solution, the mean longitudinal velocity 
U satisfies 

Proof. It follows directly from the definition of U and from the Cauchy-Schwarz and 
Holder inequalities that 



U< m \, 2 (lul 2 ) 1 / 2 . (6.2) 

L r L j h 1 / 2 ° 



Now, by (13. 3p and (13. 6p . we can estimate the term (|u| ) as follows 

Substituting (I6.3P into ( 16.21) . we obtain the result. □ 
A lower bound for the skin friction coefficient, Cf, follows directly from the theorem 
above 
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Corollary 6.2. For every stationary statistical solution, the skin friction coefficient, 
Cf, satisfies 

C >_ »g. (6.4, 

Proof. It follows immediately from Theorem 16.11 that 

Ph Ph 6iyvL x j- v 

Cf = L^P~ LJJ^WP = 2nV3 hU> 

and the result follows from the definition of the Reynolds number (14.21) . □ 
Now, we give a lower bound estimate for the mean longitudinal velocity U following 
the calculations of [I], but avoiding using an equation for the fluctuation v. 

Proposition 6.3. For every stationary statistical solution, the mean longitudinal 
velocity U satisfies 

u a sup { Tlk p ~Tk£ (5« - 2k < ft - ^ U £ " 

where 

U = {U G V; V(x, y, z) = (Utiz), 0, 0), U\ G H*(0, h), (H v (u - U)> > 0} , 

and 

II — Till 2 

H v (u - U) = J U 11 + 6(u - U, U, u - U). 

Proof. Let U G V be of the form U(x, y, z) = (Ui(z), 0, 0). We have 

(||u-U|| 2 ) = (||u|| 2 )-2((( U ,U))) + ||U|| 2 . (6.5) 
Now, since U G V is fixed, we can multiply it with the Reynolds equations and obtain 

z/(((u,U))> = (f P ,U) - (6(u,u,U)). (6.6) 
Substituting (16.61) into (16. 5p . we obtain 

^(||u|| 2 ) =u(\\u- U|| 2 ) + 2 ((f P , U) - (6(u, u, U)» - ^(||U|| 2 ). (6.7) 

Since U(x, y, z) = (Ui(z), 0, 0), by the anti-symmetry property of the trilinear term 
and by ( I4.4p . we have 



L x L y hU> -^(\\uf) 



V L X ... 2\ 

l|2 ( 6 ' 8 ) 

= + 6(u, U, u)> + 2^(f P , U) - v -h- IIUII 2 . 

Due to the form of U, we have 6(U, U, u) = 0. The orthogonality property implies 

6(u-U,U,U) =0. 
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Thus, 

6(u,U,u) = 6(u- U,U,u- U), VueV. 
Since /i is carried by V, we find 



(L x L y h) U 



= ^>^* + 6(u - U, U, u - U)> + ^(f P , U) - ^ ||U|| 2 . (6.9) 
Thus, considering only background flows U such that (H\j(v)) > 0, we have 

(L x L y h) U > ^(f P , U) - ^ ||U|| 2 . (6.10) 
We obtain the result by completing the squares. 

Remark 6.1. This theorem was shown in [I], in the context of long time averages. 
It was obtained from a derivation of an energy equation for the fluctuation variable 
v = u-U: 

\j t Mo + » ||v|| 2 + i/((v, U)) + 6(U, U, v) + 6(v, U, v) = (f P , v) (6.11) 
and by considering an energy equation for u: 

^-\u\l + v\\u\\ 2 = {L y h)PU. (6.12) 

Taking the long time average in both sides of (I6.12p . considering the same hypothesis 
for U, and substituting it into (16. lip , they have obtained the corresponding result for 
long time averages. 

However, since we want to consider any stationary statistical solutions and gen- 
eral weak solutions of the Navier-Stokes equations, we treat carefully the fluctuation 
component and avoid the energy equation (I6.1ip . The slightly modified and simpler 
proof presented in Theorem 16.31 achieves this aim. 

Theorem 6.4. For every stationary statistical solution, the mean longitudinal velocity 
and the skin friction coefficient satisfy 



2h2 - P , , f0 <p< 27 ^ VL > 



U > ( 



TJvL x ' Ah 3 ' 

2 5 / 4 vr/i 1 / 2 1/2 V2n 2 u „ 27V2n 2 u 2 L x ((> '''^ 
P 1 / 2 — if P > 

{ 3 3 / 2 j}J 2 2 W 1 Ah 3 



and 
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Cf<< 



27 1 
~2~Re : 



if < P < 



27 v / 2Vz/ 2 L a 



8tt 2 



2 Re 



Ah 3 

27v / 2vr 2 z/ 2 L x 
4^3 



(6.14) 



Proof. Following Constantin and Doering, consider the background-flow of the form 
U(x,y,z) = (E/i(z),0,0) with 



Ui(z) = { 



V 
J 

v, 

V 



0<z<5; 
5<z<h-5; 



— (h — z) , h — 5 < z < h. 
v o 



We will verify that this flow satisfies the spectral constraint H v (v) > 0, for appro- 
priate choices of V and 5. For that purpose, we bound the integral of U[(z)v 1 v 3 in 

1 1 2 

terms of 5 and ||v|| . First, divide this integral into two parts, one from to 5, and 
the other from h — 5 to h. 

In order to bound the first integral, consider the spaces H = L 2 (0,5), with the 
usual L 2 inner product, and V = {u e f/" 1 (0, 5); m(0) = 0}, with the inner product 

S ~ ~ ~ 

((u,v)) = j Q u'(z)v'(z)dz. Consider the operator A : V — > if defined by 



(Au,u) = ((u,u)), Vf G V", 



and D(A) = |u e V; iueffj. One can show that A is self-adjoint and invertible, 

with compact inverse, and that the smallest associated eigenvalue is Ai = tt 2 /AS 2 . 
Therefore, 



Ai 



\u(z)\ dz < 



r 


du(z) 


/o 


dz 



dz. 



(6.15) 



A similar statement can be made for the integral between h — S and h. Thus, the 
integral of U'{z)v\Vz can be estimated in the following way 

U ' {z)v-iVzdxdydz 
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< — 



< 



V 




J 


Jo 


V 




J J 






Viv^dxdydz 



a- 



+ 



V 



o Jo 



< 


2V5 




dv\ 




a 




7T 2 


dz 




2V5 




dvi 


< 




a 




7T 2 


dz 



* 1 
+ - 

a 
2 



\V3\ 

2a 



v i\ 

2 N 



Jo 
■dxdydz 

|2 



ViV^dxdydz 



h-8 



+ 



\V3\ 

2a 



-dxdydz 



dz 



1 

2a 



dz 



+ 



dvx 



dx 



+ 



dvo 



dy 



+ 



dvi 



dy 



+ 



dvo 



dx 



The last step above follows from the following inequality stated in [3]: 



dz 



< 



dz 



dvi 



dx 



dv 7 



dy 



+ 



dvi 



dy 



+ 



dv 2 



dx 



which is valid for divergence-free vector fields. 
Thus, choosing a = v2/2, we have 



U[(z)viV 3 dxdydz 

'o Jo 

Hence, ifu( v ) is bounded from below by 



^ ^T/X II l|2 
< — -VO V 



7T 



V 



V2, 



Therefore, Hu is non-negative if 8 < vn 2 /2y/2V, with V sufficiently large to fulfill 
the compatibility hypothesis 5 < h/2. Now, by substituting U in (I6.10p . we give a 
lower bound for U : 

vLJ/ 2 * 



U > 2L x L y [hV-SV 



5P 



(6.16) 



We maximize the lower bound above, respecting the compatibility hypotheses, with 
the following choices of V and 5: 



V 



71 



h 1 ' 2 



and 



3 1 / 2 2 3 / i L 1 x /2 



V 



P 1 ' 2 , 8 



h 2 P 



2 3/4/ l l/2pl/2 



if P > 27 



Ah 3 



5 = - 

9vL T 7 3 



if P < 27 



V2ir 2 is 2 L a 
Ah 3 



(6.17) 
(6.18) 
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The result follows immediately from the substitution of ( 16. 17ft and ( 16. 18ft into ( 16. 16ft . 

□ 

Remark 6.2. Theorem 16.41 gives a uniform upper bound estimate for the skin fric- 
tion coefficient for high Reynolds numbers. Even though this constant upper bound 
estimate is predicted by the Kolmogorov theory of homogeneous turbulence, it is 
known from experiments that corrections are necessary for turbulence in the presence 
of walls, see [3], HI [231 [28] . Actually, closure approximation theories establish the fol- 
lowing logarithmic friction law which has been confirmed by high Reynolds number 
pipe flow experiments: 

Cf ~ (j^p- (6.19) 

Thus, we conclude that while empirical arguments and experimental data predicts a 
logarithmic friction law, our rigorous mathematical bounds can only assert that 

The lower bound for Cf is of the order of the skin friction coefficient for the plane 
Poiseuille flow; see Section |HJ 

Remark 6.3. Note also that for high Reynolds number flows, the characteristic 
background velocity which leads to the estimate above is of the order of 

V ~ (h/LjWp 1 * 

while the corresponding "boundary layer" length is of the order of 

5~(vL l J 2 /h l / 2 )p- 1 / 2 . 

7. Other Estimates 

We start by deriving a lower bound for the energy dissipation rate e. 

Theorem 7.1. For every stationary statistical solution, the energy dissipation rate 
satisfies 

' 2 h 2 n2 n ^ „ ^ 27V2n 2 u 2 L x 



e > < 



27^ P ' * /0<P ^ 4.3 

2 5 / 4 vr/i 1 /2 ^ 2 v , „ 27 v / 27rVL, 

p 6 / 1 : if p > — _ _ 

I 33/2 .3/2 2 h' 3 4/i 3 



(7.1) 



Proof. The result follows from noticing that the estimate (I6.8P obtained in Proposition 



16.31 is actually a lower bound for (||u|| 2 ), and, therefore, we can follow the subsequent 
calculations in the exact same way with this term instead of U. □ 

Now, we give a lower bound on the mean kinetic energy. 
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Theorem 7.2. For every stationary statistical solution, the mean kinetic energy e 
satisfies 

e >A P _ W^ + ^V'V/^ (72) 

- 6L X L,h'\ P 12iJ + h?- 1 ' 

Proof. Taking the inner product with A~ x ip in the Reynolds equation yields 

\A- l 'H P \l = u((u, f P )> + b (u, u, A-^p) < i/<|u| ) |fp| + (\b (u, A _1 fp, u) |), (7.3) 
and since 



(|6(u,A- 1 f P ,u)|) = ( 



<2X-(/ |/»-2«||«s||tii|d2:> 



<^(/ N 2 + M 2 ^> 



(7.4) 



we find from (17. 3D that 

|^-" 2 f P |^ < f|f P | (K) 1 ' 2 + (|!> (u, A-'fp.u)!,,) 

S^dHr + f^M:) (7-5) 

Then, using ( 13. 7J , we have 

which is of the form ar 2 + br + c > for r = (||u|| 2 ) 1//2 , a = h/AL x , b = u/L x , c = 
(L y h 3 )/12L X . It gives us r 2 > b 2 /2a 2 — b{b 2 + 4ac) 1//2 /a + c/a, which implies 

2 fe 2 L 8uL x f u 2 L 2 y h 2 L^_\ 1/2 S^L.L, 



h 2 \ P 12L X J h 

□ 

Proposition 7.3. For every stationary statistical solution, the mean kinetic energy 
e satisfies 

Proof. This follows directly from inequality (16. 3p . □ 
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Theorem 7.4. For every stationary statistical solution, the energy dissipation rate e 
satisfies 



V3h 2 
6nuL 2 . 



e < --—P 2 . (71 



Proof. This follows directly from inequalities (14.41) and (16.11) . □ 
Now, we state a partial rigorous confirmation of the Kolmogorov dissipation law in 
terms of U. 

Proposition 7.5. For every stationary statistical solution, and sufficiently large pres- 
sure drop P, namely P > 27\/27i 2 u 2 L x /Ah 3 , the associated energy dissipation rate 
satisfies 

( 3 27tt 1 277r^J_\ 
' " \2^ 2 ii 2 + 4 Re + 2 7 /2 Re 2 J h' { } 

Proof. Taking the square of both sides of the second inequality in (I6.13p . we have 

3L X „ 27L x nu 27n 2 L x u 2 
2 5 /% 2 /i Ah 2 2 7 / 2 h 3 v ' 

Substituting ( 17. lOf) into (14.41) . we obtain the result. □ 

Remark 7.1. Note that we cannot invoke Proposition 15.11 neither Proposition 15.21 
to state the results (17.91) and (17. 8p in terms of their time averages. However, we can 
improve these results as follows. 

Let u(x, t) be a weak solution of the Navier-Stokes equations. It follows from the 
classical energy inequality for weak solutions of the NSE defined on [0, oo), see [5l [25] . 
that 

1 C T 1 f T 

limsup— / i/||u(s)|| ds<liminf — / (fp,u(s))ds. (7-H) 

T^oo T Jq T^oo T J 

Then, an inequality similar to (14.41) can be stated 

limsup(||u(s)|| 2 ) T < ^-liminf U T . (7.12) 



L x L y h t^oo L x T- 

Hence, if we consider inequality (16. ip for the stationary statistical solution hq that 
extends the inferior limit of the time-averaged longitudinal velocity of u(x, t), we have 



v 



^--limsup(||u( S )|| 2 ) T < f liminf U T = ^ < ^P 2 . 
L x L y h t^oo L x t^co L x buL x 

Similarly, considering inequality (16.131) for the same above, we have 

p < 3Lx TJ2 , 27 L x™ 27tt 2 L x u 2 

- 2 5 / 2 w 2 h M Ah 2 M 2 7 / 2 h 3 

3L X , a 27L x iiv n . 27n 2 L x u 2 
- (hm mt Ut) H ttt, — (hm mt U T ) 



(7.13) 



2WH 2 h K T^oo 11 Ah 2 V T^oo 11 2 7 / 2 /i 3 
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Substituting ( 17. 13ft into ( 17. 12ft . we obtain 



lim sup( 

L x L v ll T^oo 



us 



2 )t< 



+ 



27vr 1 27tt 2 1 
+ 



2 5/2 n 2 h 4 Re T 2 7 / 2 Re 2 



lim inf-r^oo Uj> 



(7.14) 



8. The plane Poiseuille flow 

In order to see how sharp our estimates are, we calculate the characteristic quan- 
tities for a specific explicit flow. It can be easily shown that there exists an explicit 
solution for the stationary version of the channel flow problem in this geometry, known 
as the plane Poiseuille flow: 

P 

upoiseuiiicO, y, z) = i^r~ z ( h ~ z ) e i- 
A straightforward calculation gives us the following estimates: 

Proposition 8.1. We have, for the plane Poiseuille flow, 

I r p2 ^4 

ep ~ lle = 2Lj^h L kw/ {h ~ zfdxdydz = 2mwf ] (8 ' 2) 

ePoiseuine = J q 4^ ~ ^fdxdydz = ^P 2 ; (8.3) 

12 

^/Poiseuille ~~ 1 (8-4) 

Ph 2 

^Poiseuille = ^2i/L ' 



and 



U Up h Ph 3 

Repoiseuille = = : 5 T ■ (8.6) 



Remark 8.1. Notice that the upper bound estimates for the mean kinetic energy, and 
mean energy dissipation obtained in the previous sections are sharp in the sense that 
they are of the same order (up to a multiplicative constant) as those just presented for 
the plane Poiseuille flow. They are sharp independently of the value of the applied 
pressure. As far as we know, these estimates were known to be sharp only when 
the applied pressure is low, since in this case the plane Poiseuille flow is globally 
asymptotically stable. 
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9. The rate of decrease of energy injection with respect to the 

scales of the flow 

In the classical theory of homogeneous turbulence, it is argued that for turbulent 
flows, the energy injection is concentrated on the large scale motions, whereas the 
energy dissipated into heat due to the molecular viscosity occurs on scales that are 
much smaller than those. 

In 1941, Kolmogorov [18] proposed that within a certain range of scales much 
lower than the energy injection scales and greater than the energy dissipative scales, 
the energy is transferred to the small scales at a nearly constant rate equal to the 
energy dissipation rate. This mechanism is called the energy cascade, and sufficient 
conditions were rigorously derived in [131 [141 [15] for the existence of this phenomenon. 
This theory was proposed in the idealized case of locally homogeneous turbulence, 
away from the boundaries, with the injection of energy restricted to the large scales. 
However, it is well known from experiments that for wall bounded turbulence, this 
hypothesis needs to be corrected; see [U[28j[23]. I n particular, the energy injection 
occurs at arbitrarily small scales. Our next result gives an upper bound for the rate 
of decrease of mean energy injection at progressively small scales. 

By taking the scalar product of the Navier-Stokes equations with the component 
u k',k" of the flow we find the energy equation for the scales of motion in the range 



^l u K ',«Ho + HI u kV"II 2 + &( u ; u > u «',«") = ((fp) «',«», iv >ls «), (9.i) 

if k" < oo, and 

2^jl u ft',oo|o + HIUk'.ooH 2 + K u 7 u > u k',oo) < ((fp)«',ooj U/e'.oo)) (9-2) 

if k" = oo. 

The energy injection at each scale associated with a wavenumber k due to the 
pressure gradient is thus given by 

±J X ±Jyll 

The energy injection in a range of wavenumbers [k', k") is given by 

i?k',k"( u ) = T — p-r((fp)/t',K") u k',k")- 
L x L y ri 

In particular, the energy injection into the wavenumbers larger than or equal to a 
given wavenumber k is given by ^ Ky00 (u). 

The mean energy injection is given by the average value of those quantities with 
respect to a given stationary statistical solution. In order to estimate the mean energy 
injection at different length scales let us prove the following lemma. 
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Lemma 9.1. The forcing term component (fp) K , for a given wavenumber k, satisfies 
l(fp)Jo = 



if k = — , k eN, k odd; 



otherwise. 
Proof. First, notice that we can write 

Aj,I,fc=K 2 

so that by the Parseval identity we have 

l(fp)Jo= E K f P.w iAfc )| 2 

We also notice that each projection satisfies 



(fp,w iM ) = f w), hk (h 
i-'x Jn 



(9.3) 



(9.4) 



(9.5) 



(9.6) 



where u> j l k denotes the first component of the eigenvector w^fc. Now, by inspecting 
the expression f!2.4[) for Wj t i )k , we notice that the integral (19.61) vanishes for all (j, I) ^ 
(0,0). Thus, ([IS]) reduces'to 



p2L * E [ h <o, k ( z ) dz 

_„2 -/0 



(9.7) 



l(fp)Jo 

Furthermore, one can deduce from the Stokes problem and the expansion (12 .4p that 
the component Wq q^z) satisfies following one-dimensional eigenvalue problem: 



515 = X 0fi,kW n k {z) 



9 2 w\ 0}k (z) 



dz 2 

(z) = 0, 



Vo,^o,o,(i( z ) 



And the normalized solution to this equation is 

1 



wo,o,ifc(0, 0, z) 



Ll /2 Ll /2 hV* 



,k7i ,k7i 
sin(— z),sm{— z),0 
n n 



with Ao,o,fc = {kn/h) 2 . Then, we can once more reduce (19.51) to 



(fp 



) I 2 



P 2 Ll 



if k = — — , for some fcG N, 
h 

if k ^ — , for every fcGM 



(9.8) 
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Then the result follows directly from the following calculation: 



Wo,o,fc(^)^ = r i/ 2r i/2, 1/2 / sm( — ^ 



2A1/2 U. if<;i8 odd, < 9 - 9 > 



0, if fc is even. □ 

We now estimate the energy injection at a given wavenumber k. 

Proposition 9.1. The mean energy injection at a given wavenumber k with respect 
to an arbitrary stationary statistical solution satisfies 

'|u K ||, it K = "i - , G JN, /e odd, 



(».(»)>< d /2 4 /2 ^« ! " "' *' ' ' (9.10) 

[0, otherwise. 

Proof. We have 

1 1 11 

(SM) = T T 7 ((fp) K ,u K ) < U~ 1/2 (f P ) K | ||u K || = -|(f P ) K | ||u K ||, 

±J X ±Jyfl J-JxlJylX LJ X LjyUj K 

and the result follows from using Lemma 19.11 □ 

Proposition 9.1. The mean energy injection on the modes larger than or equal to a 
given k and with respect to an arbitrary stationary statistical solution satisfies 

1 / 2P \ 1 / 2P \ 

(^oo(u)) < ^ ( 7ll/2Ll /2 Ll /2 h ) <Koo||> < ^ ) £l/2 " ( 9 - H ) 



Proof. We have: 



T ^ T <(fp,U*oo)> < r ^ T l^ 1/2 (fp) K ,oo|o(||u K ,oo||). 

±J X Liyll L/ X Lly!l 

Estimating the term |v4~ 1//2 (fp) Kj00 |o, we have by the Parseval identity that 

oo oo 1 



K' 

K'=K 



Then, if /c is the smallest odd number such that k < kw/h, and using Lemma I9TTI the 
following estimate holds 

AL„h 3 P 2 A 1 4L„/i 3 P 2 Z" 00 1 , 

: as 



,4-1/2,.* ,2 iLyh P ^ 1 4L^ / 



j=fe, j odd 



AL„h 3 P 2 1 ^ 4L„P 2 1 



7T 4 L X . 3/c 3 7rL;c K 3 
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Thus, 



1 2P 

-((f P ,U K ,oo)) < 



L x L y h ~ K '°°>I - ff l/2 L 3/2^/2 fc rf/2X"^~"/' 

which completes the proof. □ 
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